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Abstract. We describe families of polynomials arising in the study 
of the universal central extensions of Lie algebras introduced by Date, 
Jimbo, Kashiwara, and Miwa [DJKM83] in their work on the Landau- 
Lifshitz equations. Two of the families of polynomials we show satisfy 
certain forth order linear differential equations, are orthogonal and are 
not of classical type. 



1. Introduction 

Date, Jimbo, Kashiwara and Miwa [DJKM83] studied integrable systems 
arising from Landau-Lifshitz differential equation. The hierarchy of this 
equation is written in terms of free fermions on an elliptic curve. The au- 
thors introduced a certain infinite-dimensional Lie algebra which is a one 
dimensional central extension of q (g> C[t, t , u\u 2 = (t 2 — b 2 )(t 2 — c 2 )] where 
b 7^ ±c are complex constants and Q is a simple finite dimensional Lie alge- 
bra. This Lie algebra which we call the DJKM algebra, acts on the solutions 
of the Landau-Lifshitz equation as infinitesimal Backlund transformations. 

The Lie algebra above is an example of a Krichever-Novikov algebra (see 
([KN87b], [KN87a], [KN89]). A fair amount of interesting and fundamental 
work has be done by Krichever, Novikov, Schlichenmaier, and Sheinman 
on the representation theory of these algebras. In particular Wess-Zumino- 
Witten-Novikov theory and analogues of the Knizhnik-Zamolodchikov equa- 
tions are developed for these algebras (see the survey article [She05], and 
for example [SS99], [SS99],[She03],[Sch03a],[Sch03b], and [SS98]). 

In [CF11] the authors gave commutations relations of the universal cen- 
tral extension of the DJKM Lie algebra in terms of a basis of the algebra 
and certain polynomials. More precisely in order to pin down this central 
extension, we needed to describe four families of polynomials that appeared 
as coefficients in the commutator formulae. In this previous work we gave 
recursion relations for these polynomials and then found generating func- 
tions for them. Two of these families of polynomials are given in terms of 
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elliptic integrals and the other two families are slight variations of ultras- 
pherical polynomials. The main purpose of this note is to describe forth 
order linear differential equations satisfied by the these two elliptic families 
of polynomials (see (3.22) and (3.24)), and show that these polynomials are 
orthogonal and nonclassical (see Theorem 4.0.3 and Theorem 4.0.4). 

2. DJKM ALGEBRAS 

Let R be a commutative algebra defined over C. Consider the left R- 
module with action f(g <g) h) = fg <8> h for /, g, h G R and let K be the 
submodule generated by the elements 1® j g — f ®g— g® f '■ Then £l R = F/K 
is the module of Kahler differentials. The element / <g> g + K is traditionally 
denoted by fdg. The canonical map d : R — > Q R by df = 1 (g) / + K. The 
exact differentials are the elements of the subspace dR. The coset of fdg 
modulo dR is denoted by fdg. As C. Kassel showed the universal central 
extension of the current algebra (8> R where g is a simple finite dimensional 
Lie algebra defined over C, is the vector space g = (g <g) R) © £l R /dR with 
Lie bracket given by 

[x®f,Y®g] = [xy]®fg + (x,y)fdg, [x®f,u] =0,[u,u'] = 0, 

where x, y G g, and u, oj' G Q} R /dR and (x, y) denotes the Killing form on g. 
Consider the polynomial 

p(t) = t n + an-it™- 1 + • • • + a 

where a% G C and a n = 1. Fundamental to the description of the universal 
central extension for R = C[t, t , u\u 2 = p(t)] are the following two results: 

Theorem 2.0.1 ([Bre94], Theorem 3.4). Let R be as above. The set 

{t- 1 dt,'F I udt, . . . ,t~ n udt} 

forms a basis ofQ, R /dR (omitting t~ n u dt if oq = 0). 

Proposition 2.0.2 ([CF11], Lemma 2.0.2.). Ifu m = p(t) and R = C[t, t _1 , u\u m = 

p(t)], then in Q R /dR, one has 

(2.1) 

n-l 

((m + l)n + im)t nM ~ l udt = - ^((m + l)j + mi)a j t i+j ^ l udt mod dR 

In the Date-Jimbo-Miwa-Kashiwara setting one takes m = 2 and pit) = 
(t 2 - a 2 )(t 2 - b 2 ) = i 4 - (a 2 + b 2 )t 2 + (ab) 2 with a / ±b and neither a nor 
b is zero. We fix from here onward R = C[t, t _1 ,u \ u 2 = it 2 — a 2 ){t 2 — b 2 )]. 
As in this case ao = iab) 2 , a\ = 0, 02 = — (a 2 + b 2 ), 03 = and 04 = 1, then 
letting k = i + 3 the recursion relation in (2.1) looks like 



(6 + 2k)t k udt = -2(k - ?>){ab) 2 t k - A udt + 2fc(a 2 + b 2 )t k ~ 2 udt. 
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After a change of variables, u h-> u/ab, t i— )• t/y/ab, we may assume that 
a 2 b 2 = 1. Then the recursion relation looks like 



(2.2) (6 + 2k)t k u dt = -2(k - 3)t k ~ 4 u dt + 4kct k ~ 2 u dt, 

after setting c = (a 2 + 6 2 )/2, so that = t 4 - 2ct 2 + 1. Let P k := P k {c) 
be the polynomial in c satisfy the recursion relation 

(6 + 2k)P k (c) = AkcP k ^ 2 {c) - 2{k - 3)P fe _ 4 (c) 

for k > 0. Then set 

P(c,z) := £ P fc (c)z fc+4 = ^P fc _ 4 (c)z fc . 

fc>-4 /c>0 

so that after some straightforward rearrangement of terms we have 
= J^(6 + 2k)P k {c)z k - 4cY J kP k ^ 2 {c)z k + 2 £(fc - 3)P fc _ 4 (c)z fc 

fc>0 fc>0 fc>0 

= (-2z~ 4 + 8cz~ 2 - 6)P(c, z) + (2z~ 3 - 4cz _1 + 2z)^-P(c, z) 

+ (2z~ 4 - 8cz~ 2 )P_ 4 (c) - 4cP_ 3 (c)z- 1 - 2P_ 2 (c)z~ 2 - 4P_ 1 (c)z~ 1 . 

Hence P(c, z) must satisfy the differential equation 
(2.3) 

d_ 3z 4 -4cz 2 + l _ 2 (P_i + cP_ 3 ) z 3 + P_ 2 z 2 + (4cz 2 - 1)P_ 4 

dz (C ' Z) z 5 - 2cz 3 + z {C,Z) ~ z 5 - 2cz 3 + z 

This has integrating factor 

» {z) = exp / (i-tz'+l ~ ;) dz 

= exp(-^ ln(l - 2cz 2 + z 4 ) - ln(z)) ' 



2 zVl - 2cz 2 + z 4 ' 

2.1. Elliptic Case 1. If we take initial conditions P-3(c) = P-2(c) 
P_i(c) = and P_ 4 (c) = 1 then we arrive at a generating function 

P_ 4 (c,z) := £ F -±Ac) zk+A = Y. P ~^-^> 

k>-A k>0 

defined in terms of an elliptic integral 

4cz 2 - 1 



z\ 



P_ 4 (c,z) = z\J\-2cz 2 + z 4 y 



2 (z 4 - 2cz 2 + l) 3 / 2 



dz. 



2.2. Elliptic Case 2. If we take initial conditions P- 4 (c) = P-3(c) = 
P_i(c) = and P_2(c) = 1, we arrive at a generating function defined 
in terms of another elliptic integral: 



f>_ 2 (c, ,) = *V1 + / ( ,4_ 2 J 2 + 1)3/2 
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2.3. Gegenbauer Case 3. If we take -P-i(c) = 1, and P_2(c) = P-s(c) 
P-i{c) = and set P-\(c, z) = ^P_i in _4Z ri , then 

ra>0 

P_i(c, z) = (cz -z 3 -cz + c 2 z 3 - f; cQ^ 2 Hc)z 2n+1 ) , 



k=2 



where Qn 1 ^ 2 \c) is the n-th Gegenbauer polynomial. Hence 

P_i,_ 4 (c) = P-i,- 3 (c) = P-i,- 2 (c) = P_i, 2m (c) = 0, 
P_ 1 ,_i(c) = l, 

p , x -cQn(c) 

P-l,2n-3(c) - c a_ 1 i 

for m > and n > 2 . The Q.L~ 1/2) (c) are known to satisfy the second order 
differential equation: 

(1 " ^Q^Hc) + n(n - l)Q(-^)(c) = 

so that the P-i^ := P-i i fc(c) satisfy the second order differential equation 

d 2 d 
( fi4 " c 2 )-^P-i.2n-3 + 2c(c 2 + 1)-P_ lj2n _ 3 + {-c 2 n{n - 1) - 2)P_ lj2n „ 3 = 
ere ac 

for n > 2. 

2.4. Gegenbauer Case 4. Next we consider the initial conditions P_i(c) = 
= P_ 2 (c) = P- 4 (c) = with P_ 3 (c) = 1 and set 



P_ 3 (C, 2) = J] P- 3 ,n-4(c)^ = C 2 Z - CZ 3 - Z + CZ 3 - g Q^/ 2 ) (c 

n>0 V fc=2 

where Qn ^''(c) is the n-th Gegenbauer polynomial. Hence 

f-3,-4(c) = P-3-2(c) = P_3,_l(c) = P-l, 2m (c) = 0, 
P_3,_ 3 (C) = 1, 
p , v -Qn(c) 

cr — 1 

for m > and n > 2 and hence 

(° 2 " !)-«- P -3,2n-3 + 4C— P_ 3 ,2n-3 - (n + l)(n - 2)P_ 3 ,2n-3 = 

ore ac 

for n > 2 and P_i i2n _ 3 = cP_3,2n-3 for n > 2. 

The importance of these families of polynomials come from our previous 
work describing the universal central extension of the DJKM algebra: 



z 2n 
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Theorem 2.4.1 ([CF11]). Let g be a simple finite dimensional Lie algebra 

over the complex numbers with the Killing form (|) and define ipij(c) E 
9} R /dR by 
(2.4) 

'uii+j-2 for i + j = 1,0,-1,-2 

■P-3,i+i-2(c)(w-3 + cu-i) for i + j = 2n-l>3, n G Z, 

P_3 ii+ j_ 2 (c)(ca;_3 + /or i + j = -2n + 1 < -3, neZ, 

i P_4,|j+j|-2(c)^~4 + J P_2,|i+i|-2(c)w_ 2 /or + = 2n > 2, ii 6 Z. 



77ie universal central extension of the Date-Jimbo-Kashiwara-Miwa algebra 
is the 7j2-graded Lie algebra 

g = ?° ©g\ 

where 

g° = (g^CfM^DeCwo, g 1 = (g ® C[t, t _ V)©^-4©Cw_3eCw_ 2 eCi«;_i 
u;zf/j bracket 

[x <8 £\ y <8 t j ] = [sc, y] (8i + 5 i+i)0 i(^, y)w , 

[x ® f~V y (8) t^u] = [x, y] <8 (t l+j+2 - 2ct l+j + f +J '~ 2 ) 

+ (8 i+j -2(j + 1) - 2cj5 i+jfi + (j - l)5i+j,2) 0,yVo, 

[x (8) t i_1 u, y ® = [x, y]u (8 t i+j_1 + j(x, y)^y(c). 



3. Differential equations for Elliptic type 1 and 2 

3.1. Elliptic type 1. How we arrive at the forth order linear differential 
equation that the polynomials P-4, n satisfy stems from the approach used in 
Afken's book [AW01] for finding the second order linear differential equation 
that Legendre polynomials satisfy using only the recursion relation they 
satisfy. His technique seems to be indirectly based on ideas used in the theory 
of Grobner basis. Unfortunately the calculations and relations involved are 
rather tedious. 

From now on W6 axe going to reindcx th.6 polynomials P—^ n \ 
P- 4 (c, z) = zyfl - 2czi + * / _ ^ dz = ^ P_ 4|n(c) s" 

^ ' n=0 

= l + ,4 + ^6 + 1 (32c 2_ 5) / + i| c ( 8c 2_3 ) z 10 

5 35 1(J5 



(2048c 4 - 1248c 2 + 75) 

1155 ^ ' 



6 BEN COX, VYACHESLAV FUTORNY, AND JUAN A.TIRAO 

This means that now P_4,o(c) = 1, P-4,i(c) = P-4,2(c) = f -4.3(c) = 0. 
Besides P_4 ) o(c), the first few nonzero polynomials in c are 

P M-1 P - 4C P - 32 ° 2 ~ 5 

-T-4,4 \C) — 1, -T-4,6 — — , -T-2,8 — 

16 , - v (2048c 4 - 1248c 2 + 75) 

P_ 2 io = c(8c 2 -3), P_2i2 = -- - 

' 105 v ; ' ' 1155 

and P_4 j7l (c) satisfy the following recursion: 

(3.1) (6 + 2k)P k+i (c) = 4kcP k+2 (c) - 2(k - 3)P k (c) 

Our goal in this section is to find families of linear differential equation in c 
that these polynomials satisfy. We will call the P_ 4;n (c) DJKM polynomials. 
We start off with the generating function 



P_ 4 (c,z) = zy/l - 2cz 2 + z 4 J 



4cz 2 - 1 

■ dz 



Z 2 (z i - 2CZ 2 + l) 3 / 2 



v \ ^ 2n + 1 ^ 2n - 1 / 

\n=0 n=0 / 

where Qn^(c) is the n-Gegenbauer polynomial. These polynomials satisfy 
the second order linear ODE: 

(1 - c 2 )y" - (2A + l)cy' + n(n + 2\)y = 

where the derivative is with respect to c. Thus for A = 3/2 we get 

(3.2) (1 - c 2 )(Q( 3 / 2 ))"(c) - 4c(Q(?/ 2 ))'(c) + n(n + 3)Q(?/ 2 )(c) = 0. 

Rewrite the expansion formula for P„ 4 (c, z) to get 
(3.3) 



z -\i - 2c, 2 + z r v 2 P_ 4 ( C , Z ) = y: ^-f 1 ^ 1 - E 9 t^r 



2n-l 
6 i 



n=0 n=0 



d 2 d 

and apply the differential operator L := (1 — c 2 )— 7 — 4c— to the right 

dc A dc 

hand side to get 

L (4cQ(f/ 2 )(c)) = ((1 - c 2 )^ - 4c^) (4cg( 3 / 2 )(c)) 

= -8(n + l)Q£g J (c) - 4c(n 2 + n- 2)Q^ 2 \c) 



using the identity 

(3.4) (1 - c 2 )^ (QW(c)) = (n + 2A)cgW(c) - (n + lJQ^c). 
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After some simplification we get 
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4cQ( 3/2) (c) 2n+1 \ = ^ -8(n + l)Qgff (c) - 4c(n 2 + n - 2)Q^ /2) (c) 2n+1 

dz 



z 



cz 



(z 4 - 2cz 2 + l) 3 / 2 J z 2 {z 4 -2cz 2 + I) 3 / 2 
2 d { 1 \ cz 



V (^ 4 - 2cz 2 + I) 3 / 2 y (z 4 - 2cz 2 + l) 3 / 2 



9r 1 u > 2cz 2 • lV^' /: - 



In addition we have 



^ n 2n-l 



n=0 



2n- 1 



1,2^ 9 d 7\ ^ gg^(c) 2w _! 
4 dz 2 4 dz 4 / ^ 2n - 1 

7 n=0 



1 9 d 



1 



9 



4 dz Vz 2 (z 4 - 2CZ 2 + 1)3/2 J Az ( z 4 _ 2cz 2 + 1)3/2 

^ dz. 



4 7 Z 2 (z 4 - 2CZ 2 + 1)3/2 
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Thus the right hand side of (3.3) becomes 



Z(Z 4 - 2 C Z 2 + l) 3 / 2 J Z 2 {z*-2 C Z 2 + 1) 3 / 2 

9 d ( 1 \ cz 



dz 



cz 



9c 



dz \ (z 4 - 2CZ 2 + 1)3/2 J ( z 4 _ 2cz 2 + I) 3 / 2 



r 4 - 2cz 2 + l) 3 / 2 



dz 



1 2 d 

+ -7Z 2 — 



1 



4 dz \z 2 (z* - 2CZ 2 + 1)3/2 J 4z ( z 4 _ 2cz 2 + 1)3/2 

7 /■ 1 



+ 



4 7 Z 2 (z4 _ 2cz 2 + 1)3/2 



+ 



6cz 2 (z 2 - c) 



4z(z 4 - 2CZ 2 + 1)3/2 ^4 _ 2cz2 + 1)5 /2 

2z 2 (4z 4 - 5cz 2 + 1) 



cz 



(z 4 - 2CZ 2 + 1)3/2 z3 ( _ 2cz2 + z 4 + !) 



5/2 



9 

+ - 



4 zVz 4 - 2cz 2 + 1 



P_ 4 (c,z). 



Applying the differential operator L to the left hand side of (3.3), we get 



l(z- 1 (1-2cz 2 + z 4 )- 1/2 P-^c, 
d 2 



(l-c 2 )^(z-\l-2cz 2 +z'y l / 2 P^( 

- 2cz 2 + z A ) 
z(-4c+ (3 + 5c 2 )z 2 _4 CZ 4 



c, z 



- 4c4- ( z-\l - 2cz 2 + z 4 )" 1 / 2 P_ 4 (c, 
ac V 



(1 - 2cz 2 + z 4 ) 5 /2 



-P-4(C,Z) 



-4c + (2 + 6c 2 )z 2 - 4cz 4 d 

+ Z(l - 2CZ2 + z 4 )3/2 4lC ' Zjj 

(1-c 2 ) d 2 , 

+ Z(l - 2CZ 2 + Z 4 )V2 dc 2^-4(C,^)). 
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Hence we have 

z(-4c+ (3 + 5c 2 )z 2 -4cz 4 ) -4c+(2 + 6c 2 )z 2 -4cz 4 d 

(1 _ 2cz2 ~ z4)5/2 P-4 (C,Z)+ z(1 _ 2cz 2 + z 4)3/2 d~ c ( P ~4 (c, z) ) 

(1 - C 2 ) d 2 



s(l - 2cz 2 + z 4 )V2 rf c 2 



(P_ 4 (c,z)) 



-4cz 2 - 7 9 d 

cz 2 



4z(z 4 - 2cz 2 + l) 3 / 2 dz V (z 4 - 2cz 2 + l) 3 / 2 , 

+ 4 Z dz \z\z± - 2cz 2 + 1)3/2 ) + 4 C Wz 4 -2cz 2 + l 4(C ' Zj 

as 

- fz- 2 fz 4 - 2 C z 2 + lr 3 ^ - ~ 2{4zA - 5cz2 + 1} 

dz V {Z lCZ +L) )- Z 3(Z* - 2 C Z 2 + 1)5/2 • 

As a consequence 

z(-4c + (3 + 5c 2 )z 2 -4cz 4 ) -4c+(2 + 6c 2 )z 2 -4cz 4 d 

(1 - 2 CZ 2 + Z 4 )5/2 ^ 4(C ' Z) + z(l - 2 C Z 2 + Z 4 )3/2 & (jP " 4(C ' Z)) 

(1 - c 2 ) d 2 



s(l - 2cz 2 + z 4 ) 1 / 2 c/c 2 



(P_ 4 (c,z)) 



4cz 2 + 7 / 6z 3 (c-z 2 ) \ 

4Z(Z 4 - 2CZ 2 + 1)3/2 C ( v(z 4_ 2cz 2 + 1) 5/2j 

2 (4z 4 - 5cz 2 + 1) 

z (z 4 - 2cz 2 + if' 2 
9 1 

+ T =P-Ac, Z) 

4 z^z* - 2cz 2 + 1 V 7 
which gives us 

_9 +5cz 2_^ +4c 2^ / + 5cz 6 

= ("^ 4 - 2cz 2 + l) 2 + z 2 (-4c + (3 + 5c 2 )z 2 - 4cz 4 )^) P_ 4 (c, z) 

+ (-4c + (2 + 6c 2 )z 2 - 4cz 4 )(z 4 - 2cz 2 + 1)4 (P_ 4 (c, z)) 

dc 

+ (1 - c 2 )(z 4 - 2cz 2 + l) 2 -^ (P_ 4 (c,z)) • 
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Expanding this out in detail leads to 
-9 + 20cz 2 - (15 + 16c 2 ) z 4 + 20cz 6 

= -9 + 20cz 2 - (15 + 16c 2 ) z 4 + 20cz 6 

oo 

+ J2 ( - 9P -4,n( c ) + 20cP_ 4 ,n- 2 (c) - (l6c 2 + 6) P-4,„- 4 (c) + 20cP_ 4 ,„- 6 (c) - 9P_ 4 , n __ 8 (c) 

n=8 

- 16cP: 4jn (c) + 8(1 + 7c 2 )P: 4>n _ 2 (c) - 48c(l + c 2 )P'_^ n _ A {c) 
+ 8(1 + 7c 2 )P: 4jn „ 6 (c) - 16cP: 4>n _ 8 (c) 

+ 4(1 - c 2 ) {P'\ n {c) - \cP\ n _ 2 {c) + 2(1 + 2c 2 )P%_ 4 (c) - 4cP:' 4jn _ 6 (c) + P%_ 8 (c)) 

For the time being we will suppress the indices —4 and the variable c writing 
P-4,fc(c) as Pfc. Then we have 

(3.5) 

= -9P n + 20cP n _ 2 - (16c 2 + 6) P n _ 4 + 20cP n _ 6 - 9P n _ 8 

- 16cP; + 8(1 + 7c 2 )P^_ 2 - 48c(l + c 2 )P^_ 4 + 8(1 + 7c 2 )P^_ 6 - 16cP^_ 8 
+ 4(1 - c 2 ) {P;: - AcP;:_ 2 + 2(1 + 2c 2 )P^'_ 4 - 4cP^_ 6 + i*'_ 8 ) 

We now differentiate with respect to c the recursion relations 

(3.6) (6 + 2A;)P fc+4 = ikcP k+2 - 2(k - 3)P k 
to get 

(3.7) (6 + 2A:)P^ +4 = 4fcP fc+2 + 4fccP^ +2 - 2(k - 3)P^ 

(3.8) (6 + 2k)P% +4 = SkP' k+2 + AkcP' k \ 2 - 2(k - 3)P£. 
After setting k = n — 8 in the last equation we get 

(3.9) 

= _( 2 „ - 10)P%_ 4 + 8(n - 8)P;_ 6 + 4(n - 8)cP;'_ 6 - 2(n - ll)i£_ 8 . 

So multiplying (3.5) by (n — 11) and adding it to 2(1 — c 2 ) times the above 
gives us 

= (n - 11) ( - 9P n + 20cP n _ 2 - (16c 2 + 6) P n _ 4 + 20cP n _ 6 - 9P n _ 8 

- \QcP' n + 8(1 + 7c 2 )P^_ 2 - 48c(l + c 2 )P^_ 4 - 16ci*_ 8 

+ 4(l-c 2 ) (P:-4cP^'_ 2 )) 

+ 8 (c 2 (5n - 61) + (3n - 27)) P^_ 6 - 4(1 - c 2 )(4c 2 (n - 11) + n - 17)P"_ 4 
(3.10) 

+ 8(n-14)c(l-c 2 )P^ 6 . 
Setting k = n — 8 in (3.7) we get 
(3.11) = -(2n-10)P^_ 4 + 4(n-8)P„_ 6 + 4(n-8)cP^_ 6 -2(n-ll)P^_ 8 
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Multiplying this equation by —8c and add to the previous equation we get 
= {n - 11) ( - 9P n + 20cP n _ 2 - (16c 2 + 6) P n _ 4 - 9P n _ 8 

- 16c/* + 8(1 + 7c 2 )P' n _ 2 + 4(1 - c 2 ) (i*' - 4cP;;_ 2 ) ) 

+ 16c (-3 (c 2 + l) (n - 11) + n - 5) P' n _± - 4(1 - c 2 )(4c 2 (n - 11) + n - 17)P"_ 4 

- 12c(n - 3)P n _ 6 + 8 (c 2 (n - 29) + 3(n - 9)) P' n _ % + 8(n - 14)c(l - c 2 )P^ 6 . 

Finally if we multiply the previous equation by 2 and add it to —9 
times (3.6) we get 

= 2(n - 11) ( - 9P n + 20cP n _ 2 - 16cP^ + 8(1 + 7c 2 )P^_ 2 + 4(1 - c 2 ) {P^ - 4cP^_ 2 ) ) 

+ (6(n + 7) - 32c 2 (n - 11)) P n _ 4 + 32c (-3 (c 2 + l) (n - 11) + n - 5) P^_ 4 

- 8(1 - c 2 )(4c 2 (n - 11) + n - 17)P^_ 4 

- 60c(n - 6)P n _ 6 + 16 (c 2 (n - 29) + 3(n - 9)) P^_ 6 + 16(n - 14)c(l - c 2 )P"_ 6 . 

This is an equation without the P n _8 in it. We now get rid of the P n _e's 
in them. 

After setting k = n — 6 in (3.8) we get 
(3.12) 

= _( 2 „ - 6)i*'_ 2 + 8(n - 6)i*_ 4 + 4(n - 6)cP^„ 4 - 2(n - 9)i*'_ 6 . 

Now we multiply the previous equation by — 8(re — 14)c(l — c 2 ) and add it 
to n — 9 times the equation before it, we obtain 

= 2(n - ll)(n - 9) ( - 9P n + 20cP„-_ 2 - 16cP^ + 8(1 + 7c 2 )P' n _ 2 + 4(1 - c 2 )P^) 
+ (n - 9) (6(n + 7) - 32c 2 (n - 11)) P n _ 4 

+ 16c (c 2 - l) (n 2 - 23n + 156) P"_ 2 - 32c (c 2 (n 2 - 20n + 129) + 4n 2 - 86ra + 420) P' n _ A 

- 8 (c 2 - l) (60c 2 + n 2 - 26n + 153) P"_ 4 

- 60c(n - 9)(n - 6)P n _ 6 + 16(n - 9) (c 2 (n - 29) + 3(n - 9)) P^_ 6 
From (3.7) with k = n — 6 one has 

(3.13) 

= -(2n - 6)P,;_ 2 + 4(n - 6)P n _ 4 + 4(n - 6)cP^_ 4 - 2(n - 9)P^_ 6 . 

We multiply this equation by 8 (c 2 (n — 29) + 3(n — 9)) and add it to the 
previous equation to get 

= 2(n - ll)(n - 9)( - 9P n + 20cP n _ 2 - 16cP,; + 8(1 + 7c 2 )P' n _ 2 + 4(1 - c 2 )P^ 

- 8(2n - 6) (c 2 (n - 29) + 3(n - 9)) P' n _ 2 + 16c (c 2 - l) (n 2 - 23n + 156) P"_ 2 

- 6 (80c 2 (n - 5) - 17n 2 + 242n - 80l) P n _ 4 - 32c (l5c 2 (n - 3) + n 2 - 41n + 258) P^_ 4 

- 8 (c 2 - 1) (60c 2 + n 2 - 26n + 153) P^'_ 4 
-60c(n-9)(n-6)P n _ 6 
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From (3.6) with k = n — 6 one has 

= -(2n - 6)P n _ 2 + 4(n - 6)cP n _ 4 - 2(n - 9)P n _ 6 

We multiply this equation by — 30c(ra— 6) and add it to the previous equation 
to get 

= 2(n - ll)(n - 9) ( - 9P n - 16cP^ + 8(1 + 7c 2 )P^_ 2 + 4(1 - c 2 )P") 

+ 20c (5n 2 - 67n + 252) P n _ 2 - 8(2n - 6) (c 2 (ra - 29) + 3(n - 9)) P^_ 2 
+ 16c (c 2 - l) (n 2 - 23n + 156) P"_ 2 

+ 6 (-20c 2 (n - 4) 2 + 17n 2 - 242n + 801) P„„ 4 - 32c (l5c 2 (n - 3) + n 2 - 41n + 258) P' n _ 4 

- 8 (c 2 - l) (60c 2 + n 2 - 26n + 153) P%_ 4 

The above equation now does not have the index n — 6 in it. Now we want 
to eliminate the indices with n — 4 in them. 
From (3.8) with k = n — 4 one has 

= _( 2n - 2)i*' + 8(n - 4)i*_ 2 + 4(n - 4)cP^_ 2 - 2(n - 7)P^'_ 4 . 

We multiply this equation by — 4(c 2 — l)(60c 2 + n 2 — 26n + 153) and add it 
to n — 7 times the previous equation to get 

= 2(n - ll)(n - 9)(n - 7) T - 9P„ - 16cP^) 

+ 20c(n - 7) (5n 2 - 67n + 252) P n _ 2 

+ 6(n - 7) (-20c 2 (n - 4) 2 + 17n 2 - 242n + 80l) P„_ 4 

- 32c(n - 7) (I5c 2 (n - 3) + n 2 - 41n + 258) P^_ 4 
+ 480 (c 2 - 1) (c 2 (n - 1) - n + 9) P" 

- 32 (60c 4 (n - 4) + c 2 (-2n 3 + 45n 2 - 484n + 1749) + 15 (n 2 - 14n + 45)) P' n _ 2 

- 960c (c 2 - 1) (c 2 (n - 4) - n + 8) P^'_ 2 
From (3.7) with A; = n — 4 one has 

= _( 2n - 2)P^ + 4(n - 4)P n _ 2 + 4(n - 4)cP^_ 2 - 2(n - 7)P;_ 4 . 

We multiply this equation by — 16c(15c 2 (n — 3) + n 2 — 41n + 258) and add 
it to the previous equation to get 

= -18(n - ll)(n - 9)(n - 7)P n 

+ 480c (c 2 (n 2 - An + 3) - n 2 + 4n + 29) P' n 

+ 480 (c 2 - l) (c 2 (n - 1) - n + 9) P^' 

+ 12c (-80c 2 (n 2 - 7n + 12) + 3n 3 + 70n 2 - 1049n + 2564) P n _ 2 

- 480 (2c 4 (n 2 - 5n + 4) - 3c 2 (n 2 - 8n + 7) + n 2 - Un + 45) P' n _ 2 

- 960c (c 2 - l) (c 2 (n - 4) - n + 8) P"_ 2 

+ 6(n - 7) (-20c 2 (n - 4) 2 + 17n 2 - 242n + 80l) P n _ 4 
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From (3.6) with k = n — 4 one has 

= -2(n - l)P n + 4(n - 4)cP„_ 2 - 2(n - 7)P n _ 4 . 

We multiply this equation by 3 (-20c 2 (n - 4) 2 + 17n 2 - 242n + 80l) and 
add it to the previous equation to get 

= 120(n - 4) 2 (c 2 (n - 1) - n + 9) P n 

+ 480c (c 2 (n 2 - 4n + 3) - n 2 + 4n + 29) P^ 
+ 480 (c 2 - l) (c 2 (n - 1) - n + 9) J*' 

- 240c(n - 2) 2 (c 2 (n - 4) - n + 8) P n _ 2 

- 480 (2c 4 (n 2 - 5n + 4) - 3c 2 (n 2 - 8n + 7) + n 2 - 14n + 45) P' n _ 2 

- 960c (c 2 - l) (c 2 (n - 4) - n + 8) P"_ 2 . 
This can be rewritten as 



(3.14) 

= (n - 4) 2 (c 2 (n - 1) - n + 9) P n 

+ 4c (c 2 (n 2 - 4n + 3) - n 2 + 4n + 29) P' n 
+ 4 (c 2 - 1) (c 2 (n - 1) - n + 9) J*' 

- 2c(n - 2) 2 (c 2 (n - 4) - n + 8) P n _ 2 

- 4 (2c 4 (n 2 - 5n + 4) - 3c 2 (n 2 - 8n + 7) + n 2 - 14n + 45) P' n _ 2 

- 8c (c 2 - l) (c 2 (n - 4) - n + 8) P"_ 2 . 

We have now reduced to a differential equation with only the indices n and 
n — 2. There is a bit of a trick to reduce it down to a linear ODE with only 
the index n in it. The somewhat vague idea is to find more equations to in 
order to cancel all but terms with the index n in them. 
From (3.8) with k = n — 2 one has 

= -2(n + l)P;' +2 + 8(n - 2)P^ + 4(n - 2)cP^' - 2(n - 5)P^ 2 . 

We multiply this equation by — 4c(c 2 — l)(c 2 (n — 4) + 8 — n) and add it 
to n — 5 times the previous equation to get 

= -2(n - 5)c(n - 2) 2 (c 2 (n - 4) - n + 8) P n _ 2 

- 4(n - 5) (2c 4 (n 2 - 5n + 4) - 3c 2 (n 2 - 8n + 7) + n 2 - 14n + 45) P^_ 2 
+ (n - 5)((n - 4) 2 (c 2 (n - 1) - n + 9) P n 

+ 4c (-8c 4 (n 2 - 6n + 8) + c 2 (n 3 + 7n 2 - 105n + 177) - n 3 + n 2 + 89n - 273) P' n 

- 4 (c 2 - l) (4c 4 (n 2 - 6n + 8) + c 2 (-5n 2 + 46n - 69) + n 2 - Un + 45) P" 
+ 8(n + l)c(c 2 - l)(c 2 (n - 4) + 8 - n){n + l)P^' +2 
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From (3.7) with k = n — 2 one has 

= _2( n + l)P' n+ 2 + 4(n - 2)P n + 4(n - 2)cP' n - 2(n - 5)/*_ 2 
We multiply this equation by 

-2(2c 4 (n - 4)(n - 1) - 3c 2 (n - 7)(n - 1) + (n - 9)(n - 5)) 
and add it to the previous equation to get 

= -2(n - 5)c(n - 2) 2 (c 2 (n - 4) - n + 8) P„_ 2 

+ (-16c 4 (re 3 - 7re 2 + 14re - 8) + c 2 (re 4 + 10n 3 - 171n 2 + 416re - 256) - n 2 (n 2 - 14re + 45)) P n 

- 4c(n + 1) (4c 4 (re 2 - 6n + 8) + c 2 (-7re 2 + 60n - 93) + 3 (n 2 - 12n + 31)) P' n 

- 4 (c 2 - l) (4c 4 (n 2 - 6n + 8) + c 2 (-5n 2 + 46n - 69) + n 2 - 14n + 45) P" 
+ 4(2c 4 (n - 4)(n - 1) - 3c 2 (n - 7)(n - 1) + (n - 9)(n - 5))(n + l)P^ +2 

+ 8(n + l)c(c 2 - l)(c 2 (n - 4) + 8 - n)(n + l)P" +2 

Next we get rid of the P n _ 2 term by multiplyfing (3.6) with k = re — 2; 
= -2(n + l)P n+2 + 4(n - 2)cP n - 2(n - 5)P n _ 2 , 
by — c(n — 2) 2 (c 2 (n — 4) — n + 8) and adding it to the previous equation: 

= -n 2 (4c 4 (n 2 - 6n + 8) + c 2 (-5n 2 + 46n - 69) + n 2 - Un + 45) P n 

- Ac{n + 1) (4c 4 (n 2 - 6n + 8) + c 2 (-7n 2 + 60n - 93) + 3 (n 2 - 12n + 3l)) P' n 
-A(c 2 - l) (4c 4 (re 2 - 6re + 8) + c 2 (-5re 2 + 46re - 69) + re 2 - 14n + 45) P" 

+ 2c(re + l)(n - 2) 2 (c 2 (n - 4) - n + 8) P n+2 

+ 4(n + l)(2c 4 (n - 4)(re - 1) - 3c 2 (re - 7)(n - 1) + (n - 9)(re - 5))P^ +2 
+ 8(n + l) 2 c(c 2 - l)(c 2 (re - 4) + 8 - n)P^' +2 

Letting n \— > n — 2 in the above equation we get 

(3.15) 

= -(re - 2) 2 (4c 4 (re 2 - lOn + 24) + c 2 (-5n 2 + 66re - 181) + re 2 - 18n + 77) P n _ 2 

- 4c(re - 1) (4c 4 (n 2 - lOn + 24) + c 2 (-7n 2 + 88re - 24l) + 3 (re 2 - 16n + 59)) P' n _ 2 

- 4 (c 2 - l) (4c 4 (re 2 - lOn + 24) + c 2 (-5n 2 + 66re - 18l) + re 2 - 18n + 77) P^'_ 2 
+ 2c(re - l)(re - 4) 2 (c 2 (re - 6) - re + 10) P n 

+ 4(n - l)(2c 4 (n - 6)(re - 3) - 3c 2 (re - 9)(n - 3) + (n - ll)(n - 7))P^ 
+ 8(n - l) 2 c(c 2 - l)(c 2 (re - 6) - re + 10)P" 
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We compare this to (3.14) 
(3.16) 

= -2c(n - 2) 2 (c 2 (n - 4) - n + 8) P„_ 2 

- 4 (2c 4 (ra 2 - 5n + 4) - 3c 2 (n 2 - 8ra + 7) + n 2 - 14n + 45) P' n _ 2 

- 8c (c 2 - l) (c 2 (n - 4) - n + 8) P"_ 2 
+ (n-4) 2 (c 2 (n-l)-n + 9)P n 

+ 4c (c 2 (n 2 - 4n + 3) - n 2 + 4n + 29) P^ 
+ 4 (c 2 - 1) (c 2 (n - 1) - n + 9) j£ 

Our goal is to eliminate the P"_ 2 term, then the P^_ 2 and finally the P„- 2 
to arrive at a differential equation with just the derivatives of P n in it. We 
first have to lower the degree of c in the polynomial in front of P"_ 2 - Thus 
if we multiply (3.16) by — 2c(n — 6) and add it to (3.15) we get 

= (n - ll)(n - 2) 2 (c 2 (n + 1) - n + 7) P n _ 2 
+ 4c(n - 11) (c 2 (n 2 - l) - n 2 + 33))P^_ 2 
+ 4(n - 11) (c 2 - 1) (c 2 (n + 1) - n + 7) P^_ 2 

- 8c(n - ll)(n - 4) 2 (c 2 (n - 7) - n - l) P n 

- 4(n - 11) (c 2 (n 2 - 8n + 39) - n 2 + 8n - 7) P' n 

- 32c(n - 11) (c 2 - l) P" 

\in^\\ we get 

(3.17) = (n-2) 2 (c 2 (n + l) -n + 7)P n _ 2 

+ 4c (c 2 (n 2 - 1) - n 2 + 33))P^_ 2 
+ 4(c 2 -l) (c 2 (n + l)-n + 7)P^_ 2 

- 8c(n - 4) 2 P n 

- 4 (c 2 (n 2 - 8n + 39) - n 2 + 8n - 7) P^ 
-32c(c 2 -l)P,;'. 

We lower the degree of c in the polynomial in front of P^'_ 2 another time by 
multiplying this last equation by 2c(n — 4) and add it to n + 1 times (3.14) 
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to get 

= 8c(n-9)(n-2) 2 P n _ 2 

+ 4(n - 9) (c 2 (n 2 -An + 27) - n 2 + An + 5) P' n _ 2 

+ 32(n - 9)c (c 2 - 1) P^_ 2 

+ (n - 9)(n - 4) 2 (c 2 (n - 7) - n - l) P n 

- Ac(n - 9) (c 2 (n 2 - 12n + 35) - n 2 + Yin - 3) P' n 
+ 4(n - 9) (c 2 - 1) (c 2 (n - 7) - n - l) P". 

So that if n ^ 9, 11 one has 

(3.18) = 8c(n - 2) 2 P„„ 2 

+ 4 (c 2 (re 2 - 4re + 27) - n 2 + An + 5) P^_ 2 

+ 32c ( C 2 - 1) p;:_ 2 

+ (n-4) 2 (c 2 (n-7)-n-l)P n 

- Ac (c 2 (n 2 - 12n + 35) - n 2 + 12n - 3) P^ 

+ 4 (c 2 - 1) (c 2 (n - 7) - n - l) P" 

We lower the degree of c in the polynomial in front of P^'_ 2 another time 
by multiplying this last equation by c(n + 1) and add it to —8 times (3.17) 
to get 

= 8(n-7)(n-2) 2 P n _ 2 

+ 4c(n - 7) (c 2 (n 2 - An - 5) - n 2 + An + 37) P^_ 2 

+ 32(c 2 -l)(n-7)P^_ 2 

+ c(n - 7)(n - 4) 2 (c 2 (n + 1) - n - 9) P n 

- A(n - 7) (c 4 (n 2 - An - 5) - c 2 (n 2 + 4n - 45) + 8(n - 1)) P' n 
+ 4c(c 2 - l)(n - 7) (c 2 (n + 1) - n - 9) P^' 
Which if n ^ 7, 9, 11, then we have 

(3.19) = 8(n - 2) 2 P n _ 2 

+ 4c (c 2 (n 2 - 4n - 5) - n 2 + 4n + 37) P^_ 2 

+ 32(c 2 - i)p;:_ 2 

+ c(n - 4) 2 (c 2 (n + 1) - n - 9) P n 

- 4 (c 4 (n 2 - 4n - 5) - c 2 (n 2 + An - 45) + 8(n - 1)) P^ 
+ 4c(c 2 - 1) (c 2 (n + 1) - n - 9) P" 
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We now want to get rid of the term with P"_ 2 in it. This is done by 
multiplying the previous equation by c and add it to —1 times (3.18) we get 

= 4(c 2 -l) 2 (n-5)(n + l)P^ 2 

+ (c 2 -l) 2 (n-4) 2 (n+l)P n 

-Ac(c 2 - l) 2 (n- $)(n + l)P' n 
+ 4 ( c 2_i)3 (n + i )jF y' 

Thus as c 7^ ±1 and we are assuming n 7^ — 1,7 , 9, 11 then we have 
(3.20) = 4(n - h)P' n _ 2 + (n - 4) 2 P n - Ac(n - h)P' n + 4(c 2 - 1)P%. 
If we differentiate this with respect to c we get 
(3.21) 

= 4(n - 5)P£_ 2 + (n - 4) 2 P^ - 4(n - 5)P^ - 4c(n - 5)i*' + 8eP£ + 4(c 2 - 
= 4(n - 5)P^_ 2 + (n - 6) 2 P^ - 4c(n - 7)P£ + 4(c 2 - 

Now we work on the coefficient in front of P^_ 2 to eliminate it. 

We can multiply the previous equation by — 8(c 2 — 1) and add to n — 5 
times (3.19) to give us 

= 8(n-5)(n-2) 2 P„_ 2 

+ 4c(n - 5)(c 2 (n 2 - 4n - 5) - n 2 + An + 37)P^_ 2 
+ c(n - 5)(n - 4) 2 (c 2 (n + 1) - n - 9)P n 

+ 4(c 4 (-(n - 5) 2 )(n + 1) + c 2 (n 3 - 3n 2 - 41n + 153) - 6n 2 + 24n + 32)P^ 
+ 4c(c 2 - l)(8(n - 7) + (n - 5)(c 2 (l + n) - 9 - n))P^' 
- 32(c 2 - 1) 2 P; / 

Now we multiply (3.20) by c(37 + An - n 2 + c 2 (n - 5)(n + 1)) and add it 
to —1 times the equation above to give us 

= -8(n - 5)(n - 2) 2 P„_ 2 + 8c(n - 4) 2 (n - 1)P„ 

- 8(c 2 - l)(3n 2 - 12n - 16)J* + 192c(c 2 - 1)P" + 32(c 2 - 1) 2 P"' 

Thus if c 7^ ±1, n 7^ —1, 7, 9, 11 we have 

= -(n - 5)(n - 2) 2 P n _ 2 + c(n - 4) 2 (n - l)P n 

- (c 2 - l)(3n 2 - 12n - 16)i* + 24c(c 2 - + 4(c 2 - 1) 2 P,;" 
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If we differentiate this with respect to c we get 

-2 



= -( n -5)(n-2) 2 P^ 



+ (n-4) 2 (ra-l)P n 

+ c(n 3 - 15n 2 + 48n + 16)i^ 

+ (-c 2 (3n 2 - 12n - 88) + 3ra 2 - 12n - 40) f£' 

+ 40c(c 2 - 

+ 4(c 2 - l) 2 P,f } 

Now we multiply this by 4 and add it to (n — 2) 2 times (3.20) to get 
the 4th order linear differential equation satisfied by the polynomials P n : 



(3.22) 

16(c 2 - 1) 2 P^ } + 160c(c 2 - - 8(c 2 (n 2 -An- 46) - n 2 + 4n + 22)P^' 

- 24c(n 2 -An- Q)P' n + (n - 4) 2 n 2 P n = 0. 



3.2. Elliptic Case 2. From now on we are going to reindex the polynomials 

P-2,n- 

P_ 2 (c, z) = zVl-2cz* + z* / - _ = £ P_ a , B (c>». 

^ ' n=0 

This means that now P- 2 ,2(c) = 1, P- 2 ,3(c) = P-2,i(c) = P-2,o(c) = and 
P_2, n (c) satisfy the following recursion: 

(3.23) (6 + 2k)P k+A (c) = AkcP k+2 (c) - 2(k - 3)P fe (c) 

The first few nonzero polynomials in c are 

P_ 2 , 2 (c) = 1, P_2, 6 = 1/5, P_ 2j8 = 8c/35, 

^-2,10, = (-7 + 32c 2 )/105, P_ 2 ,i2 = 8c(-29 + 64c 2 )/1155. 

so that 

P_ 2 ( C , ,) = , 2 + ^ + ^ + + 8C(64C2 ' 29) , 12 + O(Z^). 

' v ' 5 35 105 1155 v ' 

After a very similar lengthy analysis as in the previous section we arrive at 
the following result: The forth order linear differential equation satisfied by 
the polynomial P_2 n = P n is 



(3.24) 

16(c 2 - 1) 2 P^ } + 160c(c 2 - 1)P"' - 8(c 2 (n 2 -An- 42) - n 2 + An + 18)P" 
- 24c(n 2 -An- 2)P^ + (n - 6)(n - 2) 2 (n + 2)P n = 0. 
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4. Orthogonality of the DJKM polynomials 

Now we establish the orthogonality of the family of DJKM polynomials 
P-4,n( c ) an d P-2,n( c )- After shifting the indices back by 4 we obtain for 
both families 

(4.1) 2kcP k „ 2 {c) = (3 + k)P k (c) + (k- 3)P fe _ 4 (c). 

Note that all odd polynomials are zero. Set k = 2(n + 1) and q s := i-^s- 
Then we have 

(4.2) 4(n + l)cq n = (2n + 5)q n+1 + (2n - l)q n -i, 

where q s = q s (c). For q s (c) = P_ 2 ,2s(c) 

q-l = P-2,-2 = 1) 90 = P-2,0 = 0, g-l = P_2,2 = 1/5, 

32c 2 - 7 



92 = P-2,4 = 8c/35, q 3 = P_ 



2,6 



105 

For g s (c) = P_4,2s(c) one has 



(4.3) = P_4,-4 = 1, 9-1 = P-4,-2 = 0, 50 = P-4,0 = 1, 

32c 2 — 5 16 

qi = P-4,2 = 4c/5, q 2 = P-4,4 = ^ , 93 = P-4,6 = ^ C ( 8<;2 ~ 3 ) 



If q s (c) = P_2,2s(c) we want polynomials with index n giving the degree 
of the polynomial, then we will set 

q n := q n+1 , n > -1 

while we ignore the "first" two polynomials; q-\ = 1 and qo = 0. Then 



tAA\ - Mr - o /or: - 32c 2 - 7 _ 8c(64c 2 - 29) 

4.4 90 = 1/5, 9i=8c/35, q 2 = — —— , g 3 = 

1(J5 1155 



and (4.2), becomes 

(4.5) 4(n + 2)cg n+ i = (2n + 7)g n+2 + (2n + l)9n 

or 



(4.6) 4(n + 2)cq n = (2n + 7)q n+ i + (2n + 



n— 1 , 



where the polynomials g n are of degree nine and satisfy the above recursion 
relation. 

Let w = dw(x) be a weight measure on the real line. By this we mean a 
nontrivial positive Borel measure w with finite moments of all orders. One 
introduces the following inner product in the complex linear space C[x], 

ijPiQ) = / p{x)q(x) dw(x). 

A sequence of orthogonal polynomials {9 n } n >o is a sequence in C[x] such 
that degg n = n and that (q m , q n ) = for all m ^ n. 
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We recall the following result by Favard which can be found in [Fav35] or 
in Theorem 4.4 of [Ch78]: 

Theorem 4.0.1. Let {(fin^n > 0} 6e a sequence of polynomials (monic) 
where p n is a polynomial of degree n, satisfying the following recursion 

(4.7) 4> n = (x - fi n )<f> n -i - A n n _ 2 , n = 1,2,3,... 

for some complex numbers ^ n ,X n , n = 1,2,3,... and <j>—\ = 0, cfto = 1. 
Then {4> n ,n > 0} is a sequence of orthogonal polynomials with respect to a 
(unique) weight measure if and only if \i n £l and A n +i > for all n > 1. 

We have the following equivalent version of the theorem above. 

Theorem 4.0.2. Let {p n ,n > 0} be a sequence of polynomials where p n is 
a polynomial of degree n, satisfying the following recursion 

(4.8) Xp n = a n+ ip n+ i + b n p n + C n -ip n -i, 91 = 0,1,2,... 

for some complex numbers a n +i, b n , c n -\, n = 0, 1, 2, . . . and p-\ = 0. Then 
{Pn, n > 0} zs a sequence of orthonormal polynomials with respect to some 
(unique) weight Tncasure if QTid OTtly if b n G M and c n = a n+ i / for all 
n>0. 

Proof. The equivalence of this statement with the original result of Favard 
is established in the following way: Let k n = \\4> n \\ and put p n = k~^(\> n . 
Then from (4.7) we get 

x Pn = k~ k n+ lPn+l + Vn+lPn + k~ \ n+ lk n -ip n -l. 

By putting a n+ i = k~ l k n+ i, b n = fj, n+1 and c n _i = k^Xn+iK-i we obtain 
(4.8) with b n E M and c n = a n+ \ ^ for all n > 0, because 

C n = {xPn+l,Pn) = (p n +U x Pn) = a n +l- 

Conversely, let d n be the leading coefficient of p n and put <j> n = d~ l p n . 
Then from (4.8) we get 

X(j> n -l = d~_ x a n d n (j) n + 6 n _i</> n _i + d'^Cn-zdn-z^n-^- 

Since the c/> n 's are monic we have d~_]d n d n = 1. Therefore we obtain (4.7) 
with n n = b n -i G E and X n+ i = d~ x c n -\d n _x = a n c n -i = a n a n > for all 
n > 1. □ 

2n + 7 

We apply Theorem 4.0.2 to our setting. If we set a n +i = — r, b n = 0, 

4(n + 2) 

2n + 1 

c„,_i = — p r, then the hypothesis of Favard's Theorem is not satisfied. 

4(n + 2) 

But with a modification to our recursion formula we show that the q n are 
orthogonal. 

Theorem 4.0.3. The D J KM polynomials -P-2,n( c ) are orthogonal with re- 
spect to some weight function. 
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Proof. It is sufficient to check that there exist a family of orthonormal poly- 
nomials f n and constants A n such that q n = X n f n for all n 
We have the following recursion for f n : 

4(n + 2)cf n = (2n + 7)X n+1 f n+1 + (2n + l)A n _i/ n _i 



or 



for n > 1. 
Set 



_ (2n + 7)A n+ i (2n + l)A w _i J 

4(ra + 2)A„, 4(ra + 2)A„ 



(2n + 7)A ra+1 (2n + l)A ra _i 

An+1 ~ 4(n + 2)A n ' Cn -'- 4(n + 2)A n 



rpi i ~ , , - f (2n + 5)A n (2n + l)A»_ 1 

Then A n = C n _i if and only if — — — — = — - — — — — or 

4(n + l)A n _i 4(n + 2)A n 

2 = (n + l)(2n + l) 2 
n (n + 2)(2n + 5) 

Taking Ao = 1 we can find a family of constants A n satisfying this relation. 
Hence, by Theorem 4.0.2 the polynomials f n form an orthonormal family 
with respect to some measure, and therefore the DJKM polynomials are 
orthogonal. 

□ 

Given a sequence of orthogonal polynomials {q n }n>o we consider the com- 
plex linear space of all differential operators with complex coefficients on the 
real line that have the polynomials q n as their eigenfunctions. Thus 

V(w) = {D : Dq n = j n (D)q n , 7n (D) G C for all n > 0}. 

Some properties of T>(w), see [GT07]: 

(1) The definition of T>(w) depends only on the weight function uu = w(x) 
and not on the orthogonal sequence {g n }n>o- 

(2) If D € T>(w), then 



i=0 v 7 



where fi(x) is a polynomial and deg/j < i. 
To ease the notation if v £ C let 

Mi = 1/(1/ - 1) • ■ ■ (i/ - i + 1), Mo = i- 

(3) If {gn}n>o is a sequence of orthogonal polynomials and D £ V(w) is 
of the form D = £* =0 with = £} =0 /jP)^', then 

s 

5>],/f(I?). 

i=0 
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Let us consider the Weyl algebra D = {D = J2t=o M x ) {4cY : /* e C t x ]}' 
and the subalgebra 



(4) If D G T> satisfies the symmetry condition (Dp, q) = (p, Dq) for all 
p,qeC[x], then D G V(w). 

(5) For any D G T>(w) there is a unique D* G T>(w) such that (Dp,q) = 
(p,D*q) for all p, q G C[x]. We refer to D* as the adjoint of D. The map 
D i — y D* is a *-operation in the algebra T>(w), and the orders of D and D* 
coincide. 

(6) The set S(w) = {D £ D(w) : -D = Z)*} of all symmetric differential 
operators is a real form of T>(w): 



(7) -D G D(w) is symmetric if and only if j n (D) is real for all n > 0. 
In the literature a weight w is called classical if there exists a second order 
symmetric differential operator D such that Dq n = ^ n q n for all n > 0. 

Theorem 4.0.4. T/ie sequences of orthogonal polynomials (4-4) {Qn}n>0 
and (4-3) {q n }n>o Q> T £ not classical. 

Proof. In the polynomials above we replace c by x as it is more natural to 
use this later letter as our variable. 

From (6) above, to prove the theorem is equivalent to prove that in D(w) 
there is no differential operator of order 2. Suppose D G T>(w) is of order 2. 
Then we can assume that D is of the form 



for some a,b,c,e,f G C. For (4.4) the first five terms of the orthogonal 
sequence {q n }n>o are: 




V(w) =S(w)@iS(w). 



(4.9) 




1 

qo = t, 



8,r 



32x 2 - 7 



8x(64x 2 - 29) 
1155 



9_1 = 35' 



92 = 



105 



Q3 = 



(160 x 64)x 4 - (32 x 222)x 2 + (77 x 7) 
13 x 1155 



From the definition of T>{w) and (3) we know that 
(4.10) 

/ d \ ^ d 
(ax 2 + bx + c) I — I q n + (ex + f)-j-q n = (an(n -l) + bn + c + en + f)q, 




for all n > 0. 
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If we set n = we get c + / = 0. If we put n = 1 in the above equation 
we get 

(ex + /)| = (6 + c + e + /)^ 

which implies / = and so c = (as c + / = 0) and e = b + e. Hence 6 = 0. 
Then for n = 2 we obtain 

N 32x 2 -7 2 / d \ 2 ( '32x 2 - 7\ d{32x 2 -7 
(2« + 2e)— -— = ax 2 -) l—-—)+ex- 



105 \dx J V 105 / rfx V 105 

2 64 64x 

= ax h ex 

105 105 

and hence a + e = 0. Using all this information, the above equation (4.10) 
for n = 3 is equivalent to 

8x(64x 2 - 29) 2 f d \ 2 /8x(64x 2 - 29) \ d /8x(64x 2 - 



v ; 1155 \dx J V II 55 / dx \ 



1155 



8 • 64 • 6 • ax 3 8 • 64 • 3 • ex 3 8 • 29ex 
+ 



1155 1155 1155 

which reduces to 

2a + e = 0. 

Therefore o = e = which implies that D = 0. In other words we have 
proved that the only differential operators in V(w) of degree less or equal 
to two are the constants. 
For (4.3) 

A . 32x 2 - 5 16 , 2 N 

9o = l, q 1 = 4x/5, q 2 = — , q 3 = —x (8x -3). 

35 105 ' 

A similar analysis yields also D = 0. □ 
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